Theoretical and experimental study of ion flux formation in an asymmetric 
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Parameters of a high frequency capacitive discharge in argon in axially symmetric chambers of 
different geometries are studied in experiments and by means of two-dimensional kinetic modeling 
by the ParticIe-in-Cell method. It is demonstrated that a change in the ratio of the areas of the 
driven and grounded electrodes can substantially increase the ion energy on the electrode practically 
without disturbing the plasma parameters. Particular attention is paid to studying the self-bias 
voltage and the ion distribution function on the electrode for gas pressures ranging from 15 to 70 
mTorr. The results of self-consistent calculations are in good agreement with experimental data. 
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I. INTRODUCTION 

Ion fluxes in sheaths adjacent to the electrodes gen- 
erated by a high-frequency discharge are widely used in 
various plasma technologies. Ion bombardment results 
in anisotropy of etching processes and affects the growth 
rate and the structure of films grown in the discharge. 
The applied voltage drop in the sheaths is substantial, 
and a considerable part of energy imparted to the dis- 
charge can be spent on accelerating ions. 

Processes affecting the sheath formation were studied 
in detail by many researchers (see, e.g.j^). The major- 
ity of theoretical and numerical works, however, are re- 
stricted to studying a one-dimensional problem, whereas 
real discharge geometry is usually asymmetric. In this 
case, the current density and, correspondingly, the drop 
of voltage near one electrode is greater than the corre- 
sponding values in the vicinity of the opposite electrode. 
Thus, there appears a constant self-bias voltage, which 
appreciably increases the energy of ions bombarding the 
electrode with the higher current density^. 

The discharge is asymmetric because one of the elec- 
trode is connected to the grounded walls of the discharge 
chamber, and some portion of the discharge current is 
spent there. As such a discharge has a complicated two- 
dimensional geometry, the distribution of currents can be 
calculated only with two-dimensional numerical simula- 
tions. The use of the simple ratio of the electrode areas 
(sheath capacities) in analyzing the experimental data^ 
is based on the assumption that the structure and size 
of the sheaths near different surfaces are identical. This 
assumption is not substantiated. Moreover, the sheath 
thickness and the distribution of the charged particle con- 
centrations in the sheath depends on the discharge op- 
eration mode determined, among other parameters, by 
the current density^''^. The purpose of the present work 
is a numerical and experimental study of the ion flux 
formation in a low-pressure asymmetric high-frequency 



capacitive (HFC) discharge. 

The paper is arranged as follows. The setup used in the 
experiments is described in Sec. |TI1 The kinetic model 
for the two-dimensional description of the HFC discharge 
is given in Sec. IIIIl The plasma parameters in chambers 
with different geometries are compared in Sec. IIVI The 
ion flux onto the electrode is analyzed in Sec. [V] The 
calculated results for the self-bias voltage and plasma po- 
tential, and also the distribution function of the ion flux- 
for different gas pressures are discussed in Sec. I VII The 
conclusions are formulated in Sec. IVIIl 



II. EXPERIMENTAL SETUP 

The electric parameters of a 13.56-MHz capacitive dis- 
charge in argon were measured in the experiments per- 
formed in two chambers with different configurations of 
the driven electrode. The gas-discharge chambers used 
are schematically shown in Fig. [TJ Voltage is applied 
to the top electrode with an area Arf, while the bot- 
tom electrode with an area Aq is grounded. In cham- 
ber A (see Fig. [Ha)), the effective area of the driven 
electrode is substantially smaller than the grounded elec- 
trode area {SS = Arf/Ao < 1), because the side surface 
of the discharge chamber is also subjected to a zero po- 
tential. With this ratio of the electrode areas 6S, the 
discharge is visibly asymmetric at low gas pressures, and 
the potential drop in the sheath adjacent to the driven 
electrode prevails. The potential drop at the grounded 
electrode sheath is small and depends only weakly on the 
gas pressure and on the input power. This configuration 
does not allow the flux of high-energy ions from the dis- 
charge plasma to be studied, because the ion energy ana- 
lyzer is usually mounted in the high- vacuum chamber be- 
hind the grounded electrode. Such a configuration of the 
driven electrode, however, allows measurements of the 
radial distribution of the plasma concentration and also 
the electron temperature by a moving Langmuir probe. 
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FIG. 1: Cylindrical gas-discharge chambers with different ra- 
tios of the areas of the driven and grounded electrodes. Cham- 
ber A with top driven electrode and grounded bottom elec- 
trode and side wall (a) and chamber B with top driven elec- 
trode and side wall and grounded bottom electrode (b). 



In chamber B, the driven electrode is supplemented 
with a cylindrical side part, which is shielded to pre- 
vent the breakdown on the side walls of the chamber (see 
Fig. [UJb)). In this configuration, it is only the bottom 
electrode surface that is grounded, and the ratio of the 
areas of the driven and grounded electrodes is <55' > 1. 
In chamber B, in asymmetric HFC discharge with a pre- 
vailing potential drop in the grounded electrode sheath 
is formed, and it is possible to study the high-energy ion 
flux with the use of an energy analyzer. 

The gas pressure in the experiments was varied from 6 
to 70 mTorr. The radius of chamber A is R~15 cm, and 
its height is Ho—13 cm. The discharge glows between the 
electrodes with a radius of 11 cm; the distance between 
the electrodes is d—3.8 cm. The radius of chamber B 
i?=ll cm is equal to the electrode radius. 

To measure the ion flux distribution function over en- 
ergy, there is an orifice 0.1 cm in diameter at the center of 
the grounded bottom electrode; this orifice connects the 
discharge chamber with the diagnostic chamber located 
below. The diagnostic chamber is evacuated independent 
of discharge chamber evacuation, and the pressure in the 
diagnostic chamber is lower than 1 mTorr. The diag- 
nostic chamber contains a four-grid electrostatic energy 
analyzer of the confining field. We used this analyzer to 
register the energy distribution function of the ion flux 
moving from the discharge plasma toward the grounded 
electrode. 



III. KINETIC MODEL 

The system of equations in a two-dimensional model of 
an HFC discharge with cylindrical symmetry includes the 
kinetic equations for electrons and ions (which are three- 
dimensional in terms of velocity and two-dimensional in 
space) and Poisson's equation. The energy distribution 
function for electrons fe{r,v) and ions fi{r,v) are found 



from the Boltzmann equations 
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where Ve, Vi, ne, Ui, to, and M are the electron and ion 
velocities, concentrations, and masses, respectively; Je 
and Ji are the collisional integrals for electrons and ions. 

Knowing the energy distribution functions for electrons 
and ions, we can calculate the mean energy of electrons 
and ions: 



'-fe.id Ve 
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Poisson's equation describes the electric potential dis- 
tribution 



A: 
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The boundary conditions for Poisson's equation are the 
voltage U — on the grounded electrode and U — 
UQsiniujt) + Ubias on the driven electrode. The self-bias 
voltage Ubias is calculated from the condition of a zero to- 
tal current onto the grounded surfaces and surfaces with 
applied voltage. 

System (ID)-(|11) is solved self-consistently by the 
Particle-in-Cell method with sampling of collisions by 
the Monte Carlo method (PIC MCC)''. The HFC dis- 
charge operates in argon. The kinetics of electrons in- 
cludes elastic scattering of electrons on atoms, excita- 
tion of metastable states, and ionization. Emission of 
secondary electrons from the electrodes due to bombard- 
ment by ions with the secondary emission coefficient 7 is 
also considered. 



IV. COMPARISON OF DISCHARGE 
PARAMETERS IN CHAMBERS OF DIFFERENT 
GEOMETRIES 

Let us consider the plasma parameters obtained in the 
experiment and in the self-consistent numerical solution 
of system (ID)-(|11). The discharge operates in chambers of 
different geometries (see Fig. [T]) with a fixed power of 10 
W, P=5, 30, and 70 mTorr, and 7=0.1. In chamber A, 
the electrode radius is ri— 11 cm, and the chamber radius 
is i?=15 cm. In chamber B, the bottom electrode ra- 
dius is r/=10.5 cm, and the chamber radius is i?=ll cm. 
These chambers have different ratios of the areas of the 
driven and grounded electrodes SS, because the voltage is 
applied only to the top electrode in chamber A {SS < 1) 
and to the top electrode and the side walls of the cham- 
ber {SS > 1) in chamber B. The calculated distributions 
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FIG. 2: Calculated distributions of the electron concentration 
in chambers A (a) and B (b) for P=70 mTorr. 



of the electron concentration Ue in the chambers is shown 
in Fig. [2] for P — 70 mTorr. The electron concentration 
at the center of the discharge gap is 4.9 x lO^cm"'^ in 
chamber A and 4.5 x lO^cm""^ in chamber B. In chamber 
A, the concentration Ug decreases with distance from the 
axis of symmetry toward the side wall. In chamber B, 
there is a second peak of the plasma concentration near 
the electrode edge owing to enhanced ionization. 

Figure [3] shows the distribution of the electron energy 
Eg. The energy of the electrons changes from 2.9 to 3.5 




FIG. 3: Distributions of the mean electron energy for the 
same conditions as in Fig. (5] 

eV over the volume of chambers A and B. Note that 
the maximum energy of the electrons in chamber B is 
observed between the grounded bottom electrode and the 
side walls, which is under voltage. The ionization rate 
here, however, is not too high, because the gap length 



is smaller than the characteristic length of ionization by 
electrons. 

The distribution of the electric potential (f> in Fig. |4] 
demonstrates the appearance of the self-bias voltage in- 
duced by the difference in the areas of the driven and 
grounded electrodes. Figures |Ua),(c) show the distribu- 
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FIG. 4: Distributions of the electrical potential in chambers 
A (a,c) and B (b,d) for P=70 mTorr. The potential profiles 
on (c),(d) are shown for r=0. 

tion of (j) for chamber A, where the voltage is applied 
to the top electrode with a smaller area; in this case, 
the self-bias voltage is -100 V. Figures |31Jb),(d) refer to 
chamber B where the grounded bottom electrode has a 
smaller area; the potential drop here is 110 V. 

Let us compare the measured and calculated radial dis- 
tributions of the electron concentration and energy at the 
discharge gap center. Figure [SJ) shows the measured and 
calculated concentrations of electrons in chamber A for 
different gas pressures. The calculated profiles of in 
chambers A and B are plotted in Fig. \5lh). In cham- 
ber A, which has a greater radius, the plasma concentra- 
tion monotonically decreases toward the side wall. An 
increase in the electron concentration near the edge of 
the bottom electrode is observed in chamber B. As a 
whole, the change in the chamber geometry has a weak 
effect on the plasma concentration for gas pressures rang- 
ing from 15 to 70 mTorr. With increasing pressure, the 
plasma concentration increases in both chambers from 
1.8xl0^cm-3 to 4.5xl0^cm-3. 

Figure [S] shows the radial profiles of the electron tem- 
perature Te=2/3ee at the center of the discharge gap for 
different gas pressures. The mean electron temperature 
at P = 30 mTorr, which was measured near the electrode 
edge, agrees well with the numerical data. 

The numerical and experimental data in Figs. [S] and 
[S] are in good agreement and show that the discharge 
mode and its parameters remain almost the same in both 
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FIG. 5: Radial distributions of the electron concentration at 
the center of the discharge gap () in chamber A, P=5, 30, 
and 70 mTorr (symbols, experiment), P=30 and 70 mTorr 
(curves, calculation) and (b) in chambers A and B, P=15, 30, 
and 70 mTorr (calculation). The arrows indicate the edge of 
the bottom electrode. 
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The experiment in chamber B was aimed at measuring 
the maximum energy of ions on the electrode E^ax and 
the energy of ions corresponding to the maximum 
value of the ion distribution function (IDF) in terms of 
energy. As the value of the applied voltage cannot be 
accurately determined in the experiment, we varied the 
applied voltage in our calculations to obtain an IDF close 
to the experimental curve. Figure [7] shows the measured 
and calculated dimensionless IDFs at the center of the 
bottom electrode for P = 15 and 30 mTorr in chamber 
B. The measured and calculated plasma parameters for 
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FIG. 7: Ion flux distribution functions for 15 mTorr, l7o=286 
V and 30 mTorr, l7o=260 V (b) in chamber B (the experi- 
mental and numerical data are shown by crosses and circles, 
respectively). The inset shows the potential distribution at 
r=0 for P=15 mTorr (solid curve) and P=30 mTorr (dotted 
curve) . 



FIG. 6: Radial distributions of the electron temperature in 
chamber A for P=30 mTorr (symbols, experiment) and in 
chamber B for P=15, 30, and 70 mTorr (curves, calculation). 



chambers. Because of the discharge asymmetry, however, 
the potential drop in chamber B is substantially greater 
near the grounded electrode, which makes it possible to 
form a high-energy ion flux to the electrode and to control 
the distribution function of the ions. 



the variants shown in Fig. [7] are summarized in Table ID 
The input power in the experiment was 14 W for P 15 
mTorr and 19.2 W for P = 30 mTorr. 

At gas pressures of 15 and 30 mTorr, the ions expe- 
rience several collisions as they cross the sheath; there- 
fore, the ion energy corresponding to the IDF peak is 
Ep < Emax- The maximum energy of the ions obtained 
in the experiment is E„iax^^^d eV at P = 15 mTorr and 
E,nax^i9i eV at P = 30 mTorr. It is seen from Fig. [7] 
that the experimental and calculated IDFs are in good 
agreement at the voltage Uo=285 V for 15 mTorr and at 
[/o=260 V for 30 mTorr. The inset in Fig. [7] shows the 
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calculated distribution of the potential at r = 0. The 
self-bias voltage is Ubias = W8 V for both values of the gas 
pressure. The maximum value of the plasma potential is 
0p=191 V for P = 15 mTorr and 170 V for P = 30 mTorr. 
The maximum ion energy E„iax is approximately equal 
to the averaged over 10 rf cycles plasma potential with 
respect to the grounded electrode, because the character- 
istic time of the ion crossing the sheath is much greater 
than the discharge period. 

The distribution of the charge characterizing the 
sheath width near the electrodes and chamber walls is 
shown in Fig. [51 In the case considered, the sheath width 
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FIG. 8: Distribution of the positive charge for P=30 mTorr 
and (7o=260 V. 

is 1 cm for the grounded electrode and approximately 0.5 
cm for the surface with applied voltage. At the pressure 
of 30 mTorr, the mean free path of ions with respect 
to the charge exchange collisions with neutral atoms is 
approximately 2 mm; therefore, the ions participate in 
several collisions when crossing the electrode sheath. 

Let us consider the ion flux distribution along the elec- 
trode surface. Figure [S] shows the distribution of the ion 
flux along the surface of the bottom electrode and be- 
tween the electrodes at r = 0. 

The ion flux along the bottom electrode remains al- 
most unchanged at < r < 8 cm and then decreases as 
the electrode edge is approached (see Fig. Ela)). It is 
of interest to note that the ion flux onto the top elec- 
trode is little different from the ion flux onto the bottom 
electrode, though the ion energy on the top electrode is 
considerably lower. 



VI. SELF-BIAS VOLTAGE FOR VARIOUS 
DISCHARGE GEOMETRIES 

Let us consider the self-bias voltage in chambers of dif- 
ferent geometries. Figure [TU] shows the calculated plasma 
potential relative to the grounded electrode 4>p and the 
self-bias voltage Ubias as a function of the ratio of the ar- 
eas of the driven and grounded electrodes 5S. The HFC 
discharge was calculated for chamber B with different pa- 
rameters: 1) R=ll cm, d=3.8 cm, and 5S^1.72] 2) R=7 
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FIG. 9: Distributions of the ion flux along the bottom elec- 
trode and between the electrodes at r=0 (b) for 15 mTorr, 
l7o=286 V (solid curve) and 30 mTorr, C/o=260 V (dotted 
curve) . 
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FIG. 10: Plasma potential relative to the grounded electrode 
(filled symbols) and self-bias voltage (open symbols) for the 
voltage amplitude = 160 V (a), = 260 V (b), and 
l7o = 360 V (c) for P=15, 30, and 70 mTorr. 



cm, d=3.8 cm, and bS—2.\l\ 3) R=ll cm, d=7 cm, and 
(55=2.33; 4) R=6 cm, d=8 cm, and JS'=3.9. The bot- 
tom electrode and the chamber wall are separated by a 
0.5-cm gap. It is seen in Fig. [TU]that the plasma poten- 
tial increases with increasing relative area of the driven 
electrode. In turn, Uuas increases faster with 6S and 
approaches the plasma potential value at bS >4.6. A 
decrease in the gas pressure leads to an insignificant in- 
crease in the plasma potential relative to the grounded 
electrode and the self-bias voltage. 

In our calculations, we studied the dependence of the 
voltage drop in the electrode sheaths on the electrode 
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area ratio. For a symmetric discharge, the ratio of the 
areas of the driven and grounded electrodes is SS—1, and 
the ratio (j)p/ {(j)p — Ubias) is also equal to unity. Using the 
results in Fig. [TUl we obtain a scaling exponent q for 
the expression relating the voltage drop on the electrode 
sheaths to the area of the electrodes in an asymmetric 
HFC discharge 

0p/(0p - Uuas) = SS". 

Figure [11] shows the behavior of the exponent q for dif- 
ferent pressures and voltages. Note that the curve for 




FIG. 11: Scaling exponent for the ratio (j}p/{<f)p — Ubias) = ^-S*' 
versus 5S for f/o=160 V (a), C/o=260 V (b), and [/o=360 V 
(c). 

q obtained by two-dimensional kinetic calculations is a 
nonmonotonic function, which has a maximum at SS~2- 
3: qmax='L-7 for P=15 mTorr, g„oj.=1.7 for P=30 mTorr, 
and qmax varies from 1.8 to 2.1 for P= 70 mTorr. The ex- 
ponent calculated previously^ with the use of a spherical 
model was 5=2.21. The value of q varying from 1.6 to 2.1 
was also obtained in^ by means of kinetic calculations of 
an asymmetric 2-MHz capacitive discharge. 

Figure [T^] shows the energy distribution functions in 
an ion flux in chambers with different ratios of the areas 
of the driven and grounded electrodes 5S. With increas- 
ing discharge asymmetry, the plasma potential relative to 
the smaller electrode increases, which results in a greater 
maximum energy of the ions. 



VII. CONCLUSIONS 

Formation of an ion flux in a low-pressure asymmetric 
HFC discharge was studied in experiments and by means 
of kinetic numerical simulations. Two-dimensional sim- 
ulations were performed by the Particle-in-Cell method 
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FIG. 12: Energy distribution function in the ion flux on the 
bottom electrode in chamber B for different values of SS, —30 
mTorr, and ;7o=360 V. 



with collision sampling by the Monte Carlo method. In 
the experiment, the ion flux was studied by an energy 
analyzer placed behind the grounded electrode with an 
orifice in the middle. Several reactors with different ra- 
tios of the areas of the driven and grounded electrodes 
were considered to study the effect of the chamber geom- 
etry on the ion flrrx. The plasma potential relative to the 
grounded electrode and, therefore, the maximum energy 
of the ions are demonstrated to increase with increasing 
area ratio. The measured and calculated parameters of 
the plasma, such as the electron concentration and tem- 
perature, and also the ion energy distribution functions 
are in good agreement. 
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P(mTorr) Uuas (V) 


'Pp (V) 


Uo (V) 


ne(cm"^) ji 


(s-^cm-^) 


exp cal 


exp cal 


cal 


cal 


cal 


15 109 108 


189 194 


285 


2.4x10" 


2.3x10'" 


30 108 108 


191 170 


260 


3.4x10'' 


3.0x10'" 



TABLE I: Self-bias voltage Uuas, plasma potential <f)p, ampli- 
tude of applied voltage Uo, concentration of electrons at the 
center of the discharge gap rie, and ion flux onto the bottom 
electrode ji for different gas pressures. 



